Scaling of entanglement between separated blocks in spin chains at criticality 
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We compute the entanglement between separated blocks in certain spin models showing that at 
criticality this entanglement is a function of the ratio of the separation to the length of the blocks 
and can be written as a product of a power law and an exponential decay. It thereby interpolates 
between the entanglement of individual spins and blocks of spins. It captures features of correlation 
functions at criticality as well as the monogamous nature of entanglement. We exemplify invariant 
features of this entanglement to microscopic changes within the same universality class. We find this 
entanglement to be invariant with respect to simultaneous scale transformations of the separation 
and the length of the blocks. As a corollary, this study estimates the entanglement between separated 
regions of those quantum fields to which the considered spin models map at criticality. 



Correlations have long been a central object of study 
in condensed matter with attention recently drawn to en- 
tanglement - unique correlations possible only in quan- 
tum mechanics. The presence of entanglement inside con- 
densed matter systems [H, 0] is also experimentally sup- 
ported 0]. These quantum correlations become partic- 
ularly interesting at quantum phase transitions (QPT), 
which occur at zero temperature as the relative strengths 
of interactions in a many-body system are varied Q • At 
a QPT, in general, the entanglement between individual 
spins is non-zero only for very small separations between 
the spins [fj (Fig. [TJ panel (a)). On the other hand, the 
entanglement between adjacent blocks of spins (which 
cannot by definition, have a separation) diverges with the 
length of the blocks @, 0, @] (Fig. EJ panel (b)). How- 
ever, an intermediate situation where one considers the 
entanglement between blocks of spins which are separated 
(i.e., non-complementary) is an open problem. It would 
be interesting to study how this entanglement scales, i.e., 
varies with the size of the blocks A and their separation 
x, at a QPT. Here, we conduct such a study (see, Figfl] 
panel (c)) and find that both the short ranged nature 
of spin-spin entanglement and divergent nature of adja- 
cent block entanglement can be recovered qualitatively 
as limiting cases of the expression for the entanglement 
of non-complementary blocks. This can be viewed as an 
interpolation between known behaviors of entanglement 
in the aforementioned limits. Though, what we compute 
is a form of bipartite entanglement, it also captures as- 
pects of the multipartite entanglement in the system, as 
we will discuss later. 

Note that two other information theoretic quantities 
involving non-complementary blocks of a quantum many 
body system have recently drawn much attention. One of 
them, the von Neumann entropy of disjoint blocks quan- 
tifies the entanglement of these blocks with respect to 
their complement, but not that between them [Io| . 
The second, mutual information, does quantify the cor- 
relations between the blocks, but cannot be considered 
as a measure of entanglement (i.e., the "quantum" part 
of the correlations) [l,Ill|. Seperated blocks will, in gen- 
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FIG. 1: Schematic of different types of entanglement at a 
quantum critical point. Panel (a): Critical behavior of en- 
tanglement of pairs of spins, exhibiting very limited range. 
Panel (b): Entanglement entropy of two contiguous and com- 
plementary blocks of spins, obeying a universal scaling law at 
criticality. Panel (c): Entanglement between distant blocks 
of spins as a function of their size A and separation x in a 
spin chain of finite but large length N = 2A + x . 



eral, be in a mixed state. Unfortunately, this inhibits the 
use of von Neumann entropy, for which analytical meth- 
ods exist Q, to quantifiy their entanglement and more 
involved measures [H} are necessary (as acknowledged in 
Ref. 0]). Unfortunately, even for solvable free fermionic 
models there are, in general, no known analytic meth- 
ods to compute these more involved measures. Thus, we 
choose a numerical approach which is particularly suited 
to our goal (to be justified later). 

While our principal result is the functional form of the 
entanglement of non-complementary blocks, after com- 
puting the entanglement we notice as an aside, that it 
depends only on the ratio = t and therefore exhibits 
"scale invariance" . In other words, at the considered crit- 
ical points the block-block entanglement is invariant to 
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simultaneous scale transformations x — > bx, A — > 6A. 
This is commensurate with the mapping of the citical 
models to a conformal field theory (CFT) in the thermo- 
dynamic limit. This is still an interesting observation, 
because we never explicitly use results from CFT, but 
base our study on finite but large spin chains. Previous 
computations of distant block entanglement have been 
for linear chains of quantal harmonic oscillators [l3T | . En- 
tanglement between noncomplementary regions of size R 
of quantum fields have previously been found to decay 
no faster than ~ exp [-(L/R) 2 } (for Dirac fields 0) and 
~ exp[— (L/R) 3 } (for scalar fields [l5[) where L is the 
distance between the regions. 

We focus on the finite ID XY model, consisting of a 
number N of spins | arranged on a regular lattice with 
nearest neighbour interactions and subject to an external 
magnetic field. The Hamiltonian of the spin chain reads 
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Here, a a (a = x,y, z) denote the Pauli matrices and we 
assume open boundary conditions. This model Q fea- 
tures a QPT at A = A c = 1 for any anisotropy 7. The 
XX model, 7 = 0, is also critical for field values A S [0, 1]. 
The critical regions of the model can be categorized into 
different universality classes: Ising and XX universality 
class for (7 € (0, 1], A = 1) and (7 = 0, A € [0, 1]) respec- 
tively [7[. In view of these properties, we choose three 
critical points: (i) the critical Ising (7 = 1, A = 1), (ii) 
a critical XY (7 = 0.5, A = 1) and (iii) the critical XX 
at vanishing field (7 = 0, A = 0). Our results show that 
for (i)-(iii) distant entanglement of blocks is scale invari- 
ant, and further (i) and (ii) show a common universal 
behavior with respect to this entanglement measure. 

Our approach is based on the DensityMatrix Renor- 
malisation Group (DMRG) technique [Tf|. We extend 
this method to extract from the ground state wave func- 
tion \if>) the reduced density operators pse = Ti"5^-|'0)(' ! /'l 
of system (S) and environment (E) blocks, both of equal 
length A and separated by a number x of consecutive 
spins (see Fig[TJ panel (c)), as will be outlined in detail 
in the Appendix. The DMRG procedure optimises the 
decimated set of M basis states per block in such a way, 
that the maximum entanglement entropy S of all pos- 
sible bipartitions (which arise from cutting the chain at 
a particular bond I into left and right part) is retained, 
rendering the study of critical phenomena particularly 
delicate, due to the fact that here S diverges with block 
size @, 0, H[ . In Ref . [TtJ it was found that even at criti- 
cality DMRG faithfully recovers prototypical features of 
critical ground states, such as, e.g., polynomial decay of 
correlators, as long as M is chosen sufficiently large. Note 
that DMRG faithfully reproduces the actual entangle- 
ment between the two separated blocks. We are neither 
presuming a scale invariant behaviour of the entangle- 
ment a priori (l8j to compute the representation of the 
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FIG. 2: Scaling of distant block entanglement with system 
size at and in the vicinity of the critical point A — A c = of the 
transverse Ising model (7 = 1 in Eq.JTJ. The data correspond 
to a fixed ratio p = |, and different system sizes N = 32 (open 
circles), N = 64 (squares), N = 128 (diamonds) and N = 
256 (filled circles). The crossing of dashed lines highlights 
the converged, scale invariant point for N > 128, which is 
(A, A/") = (A c , 0.052) for this ratio. 



critical ground state, nor are we course-graining blocks 
of spins into effective spins Thus, the fact that en- 
tanglement is only a function of p is an outcome of our 
calculations rather than a prior input. 

The reduced density operators ps e carry the informa- 
tion on the entanglement between the blocks S and E. 
As pse is a mixed state, the block entropy is inappro- 
priate as a measure of the entanglement. We thus have 
to use the negativity Af = (YJ i |aj| — 1) where \a,i\ de- 
note the modulus of the eigenvalues of the partial trans- 
pose (pse) Ts of pse with respect to the subsystem S, 
i.e., (wf wf \p T s s E \w s k wf) = (w s k wf\p SE \wfwf) EaEg. 
(lit; 5 )} and {Iw 6 )} are the orthogonal basis states of 
S and E respectively, chosen by the DMRG procedure. 
This is a widely used genuine measure of quantum cor- 
relations (entanglement monotone [20j) and provides a 
bound to the fidelity of teleportation with a single copy 
of the state[12]. Af depends on the size of the entangled 
regions A and their separation x. Once x exceeds the 
range of pairwise entanglement of spins, an individual 
spin in S is not entangled with an individual spin in E. 
Nonzero AT (pse) will then genuinely signal multipartite 
entanglement [ll. l2l| . We find that at the critical regions 
of the XY model, Eq[TJ Af is a function of the ratio = t 
only and shows universal behaviour. We focus on those 
values of p > 0.1, to avoid corrections stemming from fi- 
nite N, which enforce a finite value on AA (in the limit of 
adjacent blocks), whereas it is expected to diverge for in- 
finite N. We also restrict to p < 3, because beyond that 
Af assumes an order of magnitude comparable to the ac- 
curacy of the numerics. We checked the convergence of 
Af with the number of kept basis states per block M, and 
are able to produce accurate data reliable at least up to 
the 5th (3rd) decimal place for calculations in models of 
the Ising (XX) universality class. 
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FIG. 3: Negativity A/" as a function of the ratio /i = ^, for the 
critical XX at vanishing field (diamonds), critical XY (open 
circles) and critical Ising (filled circles) models. The data 
for the critical XX model are fitted according to our ansatz, 
revealing the parameters h = 0.47, a — 0.96 (dashed line). 
The critical Ising model is fitted accordingly, with parameters 
h — 0.38, a = 1.68. All data correspond to selected subsets 
of possible ratios p, in favour of visibility. 



Naturally, one might ask why we use a numerical 
method when the considered models are exactly solvable 
by a transformation to free fermions, with the ground 
states being fully specified by second moments. In the 
case of their bosonic counterparts powerful methods have 
been established [l3[ to compute the negativity of sep- 
arated blocks in terms of the second moments, yielding 
numerically exact results. However, this method does 
not readily generalize to fermions. Thus with current 
knowledge, the only option would be to explicitly con- 
struct pse in the standard (spin) basis. Even after ex- 
actly knowing the second moments it is practically unfea- 
sible to write down (store for computation) pse for large 
chains. Hence, we exploit the convenient and effective 
representation {\w ) <g> of pse arising in DMRG. 

Turning to our results shown in Fig. [2 we see that for a 
fixed ratio p the negativity Af(psE) shows the distinctive 
feature of scale invariance at the critical point A = A c = 1 
of the transverse Ising model. A similar behaviour can 
be observed for the other critical models considered here 
and for other ratios p. 

Most important and interesting is to enquire how 
N(pse) depends on different ratios p in the scale free 
point (once data have properly converged with system 
size N). From FigJ3] we see that for p < 2.5 the data 
for the two models 7 = 1 and 7 = 0.5 match near per- 
fectly signalling universality. One further recognises a 
generic shape of N(pse) as a function of p for all three 
models under consideration, its decay exhibiting a poly- 
nomial onset and an exponential tail for p > 0.2. Given 
the polynomial decay of the correlation functions with 
spatial separation, it seems plausible that quantum cor- 
relations will inherit signatures thereof. However unlike 
correlation functions, entanglement is also monogamous 
[22| : the entanglement of one block with another can only 



be a fraction of its total entanglement with the rest of the 
chain. This is another way in which the bipartite entan- 
glement between blocks captures the multipartite nature 
of entanglement in a system (i.e. between the blocks in a 
three block system) . At the limit of individual spins this 
leads the concurrence to decay very fast with separation. 
Generalizing this trend, one may expect an exponential 
decay of the entanglement of blocks with separation. In 
view of these considerations, we make the ansatz 
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with real parameters h and a. This ansatz is vindicated 
by providing an excellent fit for our data, though the 
precise values for a and h inferred from the fitting are 
somewhat sensitive (in the second decimal place) to the 
chosen fitting interval. One of the intriguing open ques- 
tions is how the exponents a and h are related to known 
critical exponents. Our fit suggests, in the case of the 
XX-model, a = 0.96, h = 0.47 and for the XY mod- 
els a = 1.68, h = 0.38. Instead, if one fitted the loga- 
rithmic negativity [l^, in the present notation defined by 
Eln = log 2 (wV+l) , to the same ansatz this would lead to 
h = 0.33 ~ 1/3. The same number was also numerically 
observed [23j for finite blocks in infinite harmonic oscilla- 
tor chains. This is further confirmation of the correctness 
of our work as both models map to the same massless 
bosonic field theory in the continuum limit. Note that 
in the limit of vanishing separation to block length ra- 
tio (vanishing p) M(pse) ~ A 1 an d thus diverges as 
A h for fixed separation (the logarithmic negativity di- 
verges as h\og 2 A). Thereby qualitatively it reproduces 
the features of block entropy (similar divergent behaviour 
was reported [24| for other measures, also including multi- 
partite settings 25]), though the negativity and block en- 
tropy are not related by any known simple formula in gen- 
eral - thus one would not exactly coincide with the other 
in any limit. In the limit of large p, i.e., small blocks of 
very distant spins the exponential part will severely dom- 
inate and ensure that their entanglement is nearly zero. 
This is expected because of the limited entangling capac- 
ity of small blocks. This capacity is exhausted by being 
entangled to their close neighbors which are granted a 
larger share of the entanglement because of the nearest 
neighbor nature of interactions. 

Summarizing, we have investigated the entanglement 
between separated blocks of spins at critical points of spin 
chains. We have conjectured an ansatz for the functional 
form of the scaling of this entanglement with separation 
and size of the blocks and shown it to be an excellent fit 
for our data. This functional form involves only a ratio of 
length scales and therefore exhibits an interesting scale 
invariance, as may be expected for models which can be 
mapped to a scale free CFT. It further qualitatively en- 
compasses two known limits of entanglement of adjacent 
blocks and that of pairs of spins. We further exemplified 
invariant features of this entanglement to microscopic 
changes within the same universality class. Interesting 
open questions are relating the numerically inferred co- 
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efficients h and a to known critical exponents, and con- 
structing measurable variables defined on the separated 
blocks which will witness their entanglement. 

APPENDIX: EXTRACTING pss WITHIN DMRG 

The canonical DMRG representation of the ground 
state \ip) is given by the schematic S • • E, where the 
boldface fonts designate a decimated representation of 
system and environment block respectively, which may 
contain a large number of physical sites. The bold dots 
represent two neighboring sites, represented in the basis 
{|t,T>, IT, I), II, T), |4,l>}- The location of this pair, in- 
dexed by the connecting bond I S [1, 2, • • • , N — 1], may 
be shifted to an arbitrary position in the corresponding 
real lattice. In particular for even N we may choose I 
such that S and E represent blocks of equal number A 
of spins, i.e. the symmetric configuration. This is the 
departure point for our study of the entanglement of the 
blocks S and E, separated by a number x of lattice sites, 
i.e. initially x = 2. In DMRG the representation of S 
arises in the course of a growth procedure from the di- 
rect product of a smaller system block S' (in our case 
containing A — 1 spins) and a single site •. The prod- 
uct S' (g) • is then transformed into a decimated basis of 
dimension M, designated as (S'»). S' is again a com- 
mon description of a yet smaller block and a single site. 
This results in a nested representation S <-> ((S" •) •) 
and so forth. In a similar fashion E arises from such 
a nesting procedure, i.e. E <-> (•(•W)). As the ba- 



sis transformations are stored, we may also go back to 
a product representation, i.e. (S'») — > S' ®«. This 
allows us to successively obtain the reduced density op- 
erators of the separated blocks pse, PS'E 1 , PS"E", and 
so on at a modest computational expense: The matrices 
representing the various pse are at most of dimension 
M 2 x M 2 . The subsequent diagonalization of (pse) Ts 
limits the maximal M we can allow, i.e., in our case 
M < 60, and conversely, the maximal N that can be 
simulated for a given accuracy. A figure of merit for 
the accuracy of the DMRG procedure is the truncated 
weight e = J2i>M w i where Wj>Wi^i> j stem from 
the Schmidt decompostion \ifj) = J2iV^\ w i) ® \ w f)- 
We require e < 10 -10 (10~ 8 ) which is satisfied for chain 
lengths as large as N = 1024 (96) in the anisotropic XY 
(XX) model, implying a similar accuracy for the block 
entropy S. With this e the negativity has converged up 
to the 5th(3rd) decimal place. 

Note added. - When finalizing this paper we became 
aware of the independent work on entanglement of non- 
complementary regions of a scalar massless field [23| . 
whose results match with our results of the XX model. 
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